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a cluster algorithm to calculate fluctuations in cluster size as well as commonly 



^. used quantities like equilibrium averages and the histograms for energy and the 

order parameter. Results provide information about the variation of cluster sizes 

•i-H 

Vh ' depending on the temperature and the number of states. They also give evidence 

for first-order transition when energy and the order parameter related measurables 
are inconclusive on small size lattices. 



Keywords: Phase transition, Potts model, clusters, Monte Carlo simulation 

fThis project is partially supported by Turkish Scientific and Technical Research 
Council (TUBITAK) under the project TBAG-1141 and TBAG-1299. 



1 Introduction 

Identifying the phase structure of a statistical mechanical system has been one of the 
most challenging problems of contemporary physics. Second-order phase transitions 
are well understood in the contexts of renormalization group ||l| and finite size scaling 
theories pj. For the case of first-order phase transitions, (for a recent review, see P]) 
the difficulty arises with the systems possesing finite, but very large correlation 
length. In such systems the transition can easily be confused with second-order 
phase transition when the measurements are done on a finite system. One of the 
commonly used methods to determine the order of a phase transition is due to Challa 
et al IQ , which considers distributions of energy or order parameter and compares 
these distibutions with a single or a double gaussian for second- or first-order phase 
transition cases, respectively. This method is based on the argument that for a first- 
order phase transition the barrier between the metastable states will be widened and 
hence distinct two gaussians will be appearent with increasing lattice size. Another 
method employed by Lee and Kosterlitz |^ ^ considers the minima in free energy, 
and similar arguments are also valid for distinguishing first-order phase transition 
from a second-order one. In both of these methods the main idea is that the system 
remains in metastable states for a considerable amount of time. Even if there exist 
metastable states, the size of the system may prevent the observation of the double- 



peak behavior of the energy distribution. For the systems with finite but very large 
correlation length, where a first-order phase transition resembles second-order one, 
one needs to compare different size lattices in order to reach a conclusion on the 
order of the transition. The well-known examples of these systems are 5-state Potts 
model in 2-dimensions [^, and 3-state Potts model in 3-dimensions P, § as well 
as finite temperature SU(3) lattice gauge theory |^. 

One of the major elements of a phase transition is the formation of clusters. Observ- 
ing the clusters and their variations in size is expected to lead a better understanding 
of phase transitions occuring in the system. In the present work, the critical be- 
haviour of the two-dimensional Potts model with 2 to 10 states is studied using 
cluster algorithm. Fluctuations in cluster size as well as the histograms for energy 
and the order parameter are obtained as a function of the temperature for different 
number of states. 

Section II gives some information about the model and the algorithm. Results and 
discussions are given in section III, and the conclusions are presented in section IV. 



2 Model and the method 



The Hamiltonian of the two-dimensional Potts model [|l0[ is given by 



n = KY.^..,^r (1) 

Here K = J/kT ; where k and T are the Boltzmann constant and the temperature 
respectively, and J is the magnetic interaction between spins cij and aj, which can 
take values 1, 2, ..., q for the g-state Potts model. The order parameter {OP) can be 
defined through the relation 

OP = ''-P-^ (2) 

where p" = N°'/L^, N"' being the number of spins with a = a, L the linear size 
and D the dimensionality of the system. The average cluster size (CS) is calculated 
by taking the average over the number of clusters (Nc), 

cs = — <yci> (3) 

where Cj is the number of spins in the i*'^ cluster. The two-dimensional Potts model 
undergoes a second-order phase transition for 2, 3 and 4 states, while a first-order 
phase transition is known to occur for higher number of states]^]. Reader can refer 



to the review article by Wu |[TT[ for detailed information about the model. 



First-order transitions can be recognized by studying double-peak behaviour of the 



energy probability distribution function P{E) |Q. When the transition is first-order, 
order parameter probabihty distribution function P{OP) should also exhibit the 
same behaviour. If the energy distribution, measured on sufficiently large lattices, 
is a single gaussian at the critical point, this is a clear indication of a second- 
order phase transition, and having two gaussian distributions is an indication of a 
first-order transition as mentioned above. Even if the energy has a single gaussian 
distribution, the order parameter or the cluster size distribution may mimic a first- 
order transition. Larger lattices may lead to double gaussian shape of the energy 
distribution but measurements made on moderate size lattices may be inconclusive. 
Moreover, a double-peak behaviour in the distribution of the order parameter may 



be misleading [O . 



Another observable which can provide information about critical behavior is the 
average cluster size (CS). The study of variations in cluster size for different K and 
q should give some indications about the nature of the transitions occuring in the 
system. By these considerations, identification of the order of phase transition is 
reduced to correct identification of the various size clusters which requires an extra 
computational effort when a usual local updating algorithm is used. The reduction 
of critical slowing down by a cluster flip algorithm was flrst introduced by Swendsen 
and Wang [^] and later modifled by Wolff |T^. The cluster algorithm used in this 



work is the same as Wolff's algorithm, with the exception that, before calculating 
the observables, searching the clusters is continued until the total number of sites in 
all searched clusters is equal to or exceeds the total number of sites in the lattice. 

3 Results and discussions 

After thermalization with 10"^ — 5 x 10^ sweeps, 5 x 10^ iterations are performed for 
32 X 32 and 64 x 64 lattices of the 2 to 10-state models at different values of the 
coupling K. Majority of the runs are started using a disordered initial configuration, 
and the averages calculated for both ordered and disordered starting configurations 
are observed to be equal within statistical errors. Longer runs with up to 10^ 
iterations are done near the finite-size critical value K^ at each state q. As a result 
of these runs, equilibrium averages for energy, cluster size and the order parameter 
are obtained, specific heat and the Binder cumulant are calculated as a function of 
K for each value of q. These quantities are used as a test ground for the algorithm 
used as well as for the thermalization and the accuracy, and they show expected 
behavior for different lattice sizes and different number of states, therefore they will 
not be reported here. Kf. values, which are estimated at the peak values of the 
specific heat for the 64 x 64 lattice are close to K^ values for the infinite lattice 
within an error of about one percent. The errors calculated using jackknife analysis 



are small (less than one percent for most of the quantities, within about 2 percent 
for the specific heat) and they can be decreased by increasing the number of Monte 
Carlo iterations. 

The histograms for energy, order parameter, and the average cluster size as a function 
of the number of iterations are obtained for g = 2 to 10 near Kc- The models with 
g = 2 and 3 have similar histograms with a single gaussian, and for the cases g > 5, 
all histograms have distinct double peaks, getting higher and sharper as q increases. 
In the histograms for energy, double peaks occur for g > 5 with almost equal maxima 
near the Kc values estimated using the specific heat data, as expected. Figure 1 
shows the histograms for energy and the order parameter for g = 3,4 and 5 on 
64 X 64 lattice obtained at K^. Plots for g = 3 have single peaks for both energy and 
the order parameter. The histograms for g = 5 show indication of the first-order 
character, with the appearance of double peaks in both plots. The energy histogram 
for g = 4 is similar to the one for g = 3, but the order parameter histogram shows 
a variation which looks more similar to the plot for g = 5. 

The fluctuations PCS {PCS =< {CSy > - < CS >^) in the average cluster 
size CS are calculated for different values of g on 64 x 64 lattice as a function of 
K. Figure 2 shows PCS for g varying from 3 to 7 and the CPU times for the 
corresponding Monte Carlo runs. As can be seen from these plots, PCS and the 



CPU time give similar information about the formation of the clusters {1/FCS has 
a similar variation as the CPU time). When K is small (i.e. the temperature is 
high), the cluster sizes are small, hence FCS is small, and the algorithm spends 
a long time to search for the clusters small in size but large in number. As K 
increases, the cluster sizes are getting larger, so the CPU time decreases, but FCS 
increases due to the existence of the small-size clusters as well as the large ones. 
At the critical point, formation of the largest clusters results in the largest value of 
FCS, with the corresponding minimum value of the CPU time, since there are only 
a few clusters to be searched. One can make a similar discussion when the critical 
point is approached from above (in the region where K > Kc), with the exception 
that the CPU time is almost constant for large K, as mentioned in section II. 

FCS plots also show that the maximum value of FCS increases as q increases. 
When q is getting larger, decreasing correlation length (due to increasing first-order 
character) results in breaking large clusters into smaller ones. Near the critical point, 
the tendency of forming large clusters together with the disintegration process causes 
FCS to increase with increasing q. In the second-order phase transitions, long range 
correlations keep the clusters as they are for a very long time or for large number of 
iterations (hence FCS is small for small q). The FCS and CPU plots show that the 
models with g > 4 behave differently at two sides of the critical point. The plots for 



g = 3 are in the form of curves with a smooth variation, passing through a minimum 
near Kf.. Plots for the other cases show different behaviour in two separate regions 
around the critical point, with a sharp variation near K^. value. The CPU data for 
5 = 3 can be fitted to a polynomial function around K,, at both sides of the critical 
point, but no fit seems to be possible for g > 4. The only possibility is the fit of two 
distinct polynomials to two separate curves, for K < Kc and K > Kc- 



4 Conclusions 

The two-dimensional Potts model with 2 to 10-states is studied using cluster algo- 
rithm for lattices with sizes 32 x 32 and 64 x 64. Energy, order parameter, cluster 
size, fiuctuations in the cluster size, specific heat and the Binder cumulant are cal- 
culated, and the histograms for energy, order parameter and the cluster size are 
obtained as a function of the coupling K. 

Studies in this work are focused on the cluster size related observables. One of these 
observables is the fiuctuations in the cluster size (FCS) which measures the extent 
of the coexistence of different size clusters. This operator is very sensitive to the 
correlation length and exhibits sharp peak at the transition region if the transition 



is first-order. The similar information can also be obtained from the CPU time 
spent to calculate the averages for a certain number of iterations. If there exist only 
small clusters, due to the modified algorithm used in this work, the CPU time is 
extremely large and growing cluster size reduces this time. In the first-order case, 
small clusters grow very fast near Kc-, while in the second-order case, the average 
cluster size grows smoothly. Both figures 2. a (FCS) and 2.b (CPU) indicate that on 
64 X 64 lattice g = 3 and 4 have different character than g > 5 Potts models. From 
these two figures one can conclude that for the systems with varying correlation 
length, FCS and CPU are sensitive to the correlations even though the correlation 
length is much larger than the system size. 

Further work is planned, using the same algorithm, to study cluster distributions in 
detail as well as random-field and random-bond Potts models. 
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Figure captions 

Figure 1. a) Energy, b) the order parameter histograms, as a function of the number 
of iterations for g = 3, 4 and 5 on 64 x 64 lattice near Kc- 

Figure 2. a) Fluctuations in the average cluster size, b) the CPU time, as a function 
of K for q varying from 3 to 7 on 64 x 64 lattice. Computations are done using 
IBM-RS/6000 workstations. The errorbars in FCS plots are about the size of 
the data points and are visible only for large q. 
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